Long-range quantum entanglement in noisy cluster states 
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We describe a phase transition for long-range entanglement in a three-dimensional cluster state 
affected by noise. The partially decohered state is modeled by the thermal state of a short-range 
translation-invariant Hamiltonian. We find that the temperature at which the entanglement length 
changes from infinite to finite is nonzero. We give an upper and lower bound to this transition 
temperature. 
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I. INTRODUCTION 

Nonlocality is an essential feature of quantum mechan- 
ics, put to the test by the famous Bell inequalities Q and 
verified in a series of experiments, see e.g. 0- Entangle- 
ment |3 is an embodiment of this nonlocality which has 
become a central notion in quantum information theory. 

In realistic physical systems, decoherence represents 
a formidable but surmountable obstacle to the creation 
of entanglement among far distant particles. Devices 
such as quantum repeaters Q and fault-tolerant quan- 
tum computers are being envisioned in which the en- 
tanglement length 0, is infinite, provided the noise 
is below a critical level. Here we are interested in the 
question of whether an infinite entanglement length can 
also be found in spin chains with a short-range inter- 
action that are subjected to noise. A prerequisite for 
our investigation is the existence of systems with infi- 
nite entanglement length at zero temperature. An exam- 
ple of such behavior has been discovered by Vcrstraete, 
Martin-Delgado, and Cirac with spin-1 chains in the 
AKLT- model la, and by Pachos and Plenio with cluster 
Hamiltonians see also |lfj| . In this paper, we study 
the case of finite temperature. We present a short-range, 
translation-invariant Hamiltonian for which the entangle- 
ment length remains infinite until a critical temperature 
T c is reached. The system wc consider is a thermal cluster 
state in three dimensions. We show that the transition 
from infinite to finite entanglement length occurs in the 
interval 0.30 A < T c < 1.15 A, with A being the energy 
gap of the Hamiltonian. 

We consider a simple 3D cubic lattice C with one spin- 
1/2 particle (qubit) living at each vertex of the lattice. 
Let X u , Y u , and Z u be the Pauli operators acting on the 
spin at a vertex u G C. The model Hamiltonian is 



H 
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Here neigh (u) is a set of nearest neighbors of vertex 
u. The ground state of H obeys eigenvalue equations 
K u \4>)c = \4>)c and coincides with a cluster state [Tlj . 



We define a thermal cluster state at a temperature T as 



Pes = iexp(-/3if), 



(2) 



where Z = Tr e " H is a partition function and (3 = T 1 . 
Since all terms in H commute, one can easily get 



Pes 



L I] {I + tanh (/3A/2) K u 
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(3) 
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Let A, B C C be two distant regions on the lattice. Our 
goal is to create as much entanglement between A and B 
as possible by doing local measurements on all spins not 
belonging to AU B. Denote a as the list of all outcomes 
obtained in these measurements and B as the state of 
A and B conditioned on the outcomes a. Let E[p) be 
some measure of bipartite entanglement. Following || 
we define the localizablc entanglement between A and B 
as 



E(A, B) = nrnxJ2PaE[p-* B ], 



(4) 



where p a is a probability to observe the outcome a and 
the maximum is taken over all possible patterns of lo- 
cal measurements. To specify the entanglement measure 
E[p] it is useful to regard p* B as an encoded two-qubit 
state with the first logical qubit residing in A and the 
second in B. We choose E[p] as the maximum amount of 
two-qubit entanglement (as measured by entanglement of 
formation) contained in p. Thus < E(A, B) < 1 and an 
equality E(A, B) = 1 implies that a perfect Bell pair can 
be created between A and B. Conversely, E(A, B) = 
implies that any choice of a measurement pattern pro- 
duces a separable state. 

In this paper we consider a finite 3D cluster 

C = {u = (ui, U2, u 3 ) : 1 < m, U2 + 1 < I; 1 < "3 < d} 

and choose a pair of opposite 2D faces as A and B: 

A = {u E C : u 3 = 1}, B = {u G C : u 3 = d}, 

so that the separation between the two regions is d — 1. 
In Section HU we show that |2U| 



Ld 



lim E(A, B) = 1 for T < 0.30 A. 
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Further, we show in Section lILTl that if T > 1.15 A then 
E(A, B) = for d > 2 and arbitrarily large I. 



II. LOWER BOUND 

We relate the lower bound on the transition tempera- 
ture to quantum error correction. From Eq. @ it follows 
that p cs can be prepared from the perfect cluster state 
\4>}c by applying the Pauli operator Z u to each spin aeC 
with a probability 



1 



P = 



1 + exp(/3A) ' 



(5) 



Thus, thermal fluctuations are equivalent to independent 
local Z-errors with an error rate p. 

We use a single copy of p cs and apply a specific pat- 
tern of local measurements which creates an encoded Bell 
state among sets of particles in A and B. For encoding 
we use the planar code, which belongs to the family of 
surface codes introduced by Kitaev. The 3D cluster state 
has, as opposed to its ID counterpart [Til ], an intrinsic 
error correction capability which we use in the measure- 
ment pattern described below. Therein, the measure- 
ment outcomes are individually random but not indepen- 
dent; parity constraints exist among them. The violation 
of any of these indicates an error. Given sufficiently many 
such constraints, the measurement outcomes specify a 
syndrome from which typical errors can be reliably iden- 
tified. The optimal error correction given this syndrome 
breaks down at a certain error rate (temperature), and 
the Bell correlations can no longer be mediated. This 
temperature is a lower bound to T c , because in principle 
there may exist a more effective measurement pattern. 

To describe the measurement pattern we use, let us 
introduce two cubic sublattices T e ,T a C C with a double 
spacing. Each qubit u £ C becomes either a vertex or 
an edge in one of the sublattices T e and T a . The sets of 
vertices V(T e ) and V(T ) arc defined as 

V(T e ) = {u= (e,e,e)eC} } 
V{T ) = {u = (o, o, o) G C}, 

where e and o stand for even and odd coordinates. The 
sets of edges E(T e ) and E(T a ) are defined as 

E(T e ) = {u = (e,e,o), (e,o, e), (o, e, e) G C}, 
E{T Q ) = {u= (o,o,e),(o,e,o), (e,o,o) eC}. 

The lattices T e , T Q play an important role in the identi- 
fication of error correction on the cluster state with a Z2 
gauge model ^2| . They are displayed in Fig. [5] 

Let us assume that the lengths I and d are odd [22^ . 
The Bell pair to be created between A and B will be 
encoded into subsets of qubits 

L = {u = (o, e, 1), (e, o, 1) G C} C A, 
R = {u = (o, e, d), (e, o,d) G C} C B. 




113 even 



113 odd 



FIG. 1: (color online) (a) Measurement pattern on the first 
and last slice of C, for / = 5. The resulting state is in the 
code space of the planar code. (The unmeasured qubits are 
displayed as shaded circles.) (b) Lattice for the planar code. 
(c)/(d) Measurement pattern for even and odd inner slices. 



Each qubit u G C is measured either in the Z- or X- 
basis unless it belongs to L or R. Denoting Mx and Mz 
local X- and Z-measurements, we can now present the 
measurement pattern: 



M z : Vu G V(T e ) U V(T ), 

M x : V«e E(T e ) U E(T a )\(L U R), 



(G) 



We denote the measurement outcome ±1 at vertex u by 
z u or x u , respectively. A graphic illustration of the mea- 
surement patterns for the individual slices is given in 
Fig. H] 

Before we consider errors, let us discuss the effect of 
this measurement pattern on a perfect cluster state. Con- 
sider some fixed outcomes {x u }, {z u } of local measure- 
ments and let \iP)lr be the reduced state of the unmea- 
sured qubits L and R. We will now show that \iP)lr is, 
modulo local unitaries, an encoded Bell pair, with each 
qubit encoded by the planar code [Tflj . the planar coun- 
terpart of the toric code 0]- The initial cluster state 
obeys eigenvalue equations K u \<f>)c = \<t>)c- This implies 
for the reduced state 



Zp,uW)lr = ^p,uWilr, Vu=(e, e, 1), 



(7) 



where Z P , U = <g> t , encigh(u)nL Z v is a plaquette (z-type) 
stabilizer operator for the planar code The eigen- 

value Xp tU depends upon the measurements outcomes as 
Xp.u = ^ti x (ui,M2,2)' Note that in the planar code the 
qubits live on the edges of a lattice rather than on its ver- 
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tices. The planar code lattice is distinct from the cluster 
lattice C; see Figs. IT121 



l-l 







l-l 








FIG. 2: (color online) The measurement pattern on the clus- 
ter C. The sublattices T e and T are displayed (thick lines). 
For reference, the cluster lattice is also shown (thin lines) and 
the axis labeling shows the cluster coordinates. Cluster qubits 
measured in the Z-basis (on the sites of T and T e ) axe dis- 
played in black, and qubits measured in the X-basis (on the 
edges of T and T e ) are displayed in gray (red). The large 
circles to left and to the right denote the unmeasured qubits 
which form the encoded Bell pair. The measurement pattern 
has a bcc symmetry. 



From the equation ILeneig 
(o, o, 1), we obtain 



U(„ 



A' 



for u 



Xs,u\ 



)LR 



^S,uW)LR, Vli= (o, O, 1), 



(8) 



where X s , u = ® l)en eigh(«)nz, X v coincides with a site (re- 
type) stabilizer operator for the planar code and 
As 

n«enei g h («) z v, where neigh D refers to 
a neighborhood relation on the sublatticc T Q . The code 
stabilizer operators in Eq. Q and (jSJ) are algebraically 
independent. There are (I 2 — l)/2 code stabilizer gener- 
ators for (I 2 + l)/2 unmeasured qubits, such that there 
exists one encoded qubit on L. By direct analogy, there 
is also one encoded qubit located on R. 

Next, we show that \iji)LR is an eigenstate of X^Xr 
and ZlZr, where X and Z are the encoded Pauli op- 
erators X and Z, respectively, i.e. \iP)lr is an encoded 
Bell pair. The encoded Pauli operators ^3 on L and 
R are X L[R] = (g) X (uuU2A [ d]) for any even u 2 , and 

Ui odd 



l-l 







T nontrivial trivial 
1 o error cycle error cycle 



i "i 



d 



FIG. 3: (color online) A homologically nontrivial and a ho- 
mologically trivial error cycle on the lattice T . The nontriv- 
ial error cycle stretches from one rough face to the opposite 
one while the trivial error has both ends on the same face. 
Only the qubits belonging to T D are shown and the qubits im- 
portant for establishing the X^^s-correlation are displayed 
enlarged. 



z 



L[R] 



® ^(«i,«2,i[dl) f° r an y °dd u\. To derive the 



ii2 even 

Bell-correlations of \i(>)lr let us introduce 2D slices 



T (u 2 ) 

- 1 XX 
-L ZZ 



= {u = (o, u 2 , o) G C} C T , 
— {u = (m, e, e) G C} C T e . 



The eigenvalue equation FJ („ 2 ) A„ 
even U2 implies for the reduced state 



X L X R \4>)lr = \xx W) 



LR, 



with 



(9) 



with X X x = II 



vET 



("2 + 1) 

XX 



UT 



("2-1) ^t, 



n 



Here and thereafter it is understood that 
for all u G* C Similarly, from |0)c = J 



'\(LUR) ' 



Zu = 1 



vet: 



for u\ odd, we obtain for the reduced state 
ZlZr \4>)lr = ^zz \iP)lr, 



n 



VET. 



("i) X v . 



(10) 



Thus the 



with x zz = n„ eT (y+D UT ^i-i) «o 

eigenvalue Eqs. (|7I10|I show that the measurement pat- 
tern of Eq. JBJ projects the initial perfect cluster state 
into a state equivalent under local unitaries to the Bell 
pair, with each qubit encoded by the planar code. 

It is crucial that the measurement outcomes {z u } and 
{x v } are not completely independent. Indeed, for any 
vertex u G T with 1 < u% < d the eigenvalue equation 



n 



i;£ncigh(u) 



K v \4>)c = k 

n ^ 

i>Eneigh(u) 



)c implies the constraint 

• n Zw = l ( u ) 

u/£neigh (u) 



Analogously, for any vertex u £ T e one has a constraint 

1, (12) 



-u^ncigh(u) 



n «> 
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where neigh e refers to a neighborhood relation on the 
lattice T e . Thus there exists one syndrome bit for each 
vertex of T e and T , (with exception for the vertices of 
T a with U3 = 1 or U3 = d). 

What are the errors detected by these syndrome bits? 
Since we have only Z-crrors (for generalization, see re- 
mark 1), only the Admeasurements arc affected by them. 
Each A-mcasurcd qubit is cither on an edge of T a or 
T e . Thus, we can identify the locations of the elemen- 
tary errors with E(T a ) and E(T e ). From the equations 
Eq. (|11I12|) , each error located on an edge creates a syn- 
drome at its end vertices. 

Let us briefly compare with [l2T |. Therein, independent 
local X-and Z-errors were considered for storage whose 
correction runs completely independently. The A-errors 
in this model correspond to our Z-crrors on qubits in 
E(T e ), and the Z-storage errors to our Z-errors on qubits 
in E(T Q ), if the X— and Z— error correction phases in [l^ 
are pictured as alternating in time. 

The syndrome information provided by Eqs. i|llll2[) is 
not yet complete. There are two important issues to be 
addressed: (i) There are no syndrome bits at the ver- 
tices of T with U3 = 1 or H3 = (I; (ii) Edges of T e with 
113 = 1 or 113 = d have only one end vertex, so errors 
that occur on these edges create only one syndrome bit. 
Concerning (i), to get the missing syndrome bits we will 
measure eigenvalues Xp tU and Xs, u for the plaquette and 
the site stabilizer operators living on the faces A and 
B, see Eqs. (I7I8II . Such measurements are local oper- 
ations within A or within B, so they can not increase 
entanglement between A and B. For any u = (o, o, 1) or 
u = (o, o, d) it follows from Eq. JSJ that 

Ag, u a;( Ul ,u2,2) Zu Zv = 1, foru 3 = l, 

u£neigh (u) 

^S,uX(u u ua,d-l) z u z v = lj i0TU 3 = d. 

■u(Encigh (u) 

(13) 

For any vertex u = (0,0,1) or u = (0,0, d) there are 
several edges of the lattice T D incident to u. It is easy to 
see that a single Z-error that occurs on any of these edges 
changes a sign in Eqs. (|13[1 . Thus, these two constraints 
yield the syndrome bits living at the vertices u — (o, o, 1) 
and u = (o, o, d), so the issue (i) is addressed. Concerning 
(ii) , we make use of Eq. (Q and obtain 

Ap, u a; u Z( UliM2i2) = 1, for any u = (e,e, 1), 
Xp,u x u Z( Ul ,u 2 ,d-i) = 1, for any u = (e,e,d). 

Since we have only Z-errors, the eigenvalues Xp_ u and 
the outcomes Z( Ml . tl2 ,2)j Z(vi,u 2 ,d-i) are n °t affected by 
errors. Thus the syndrome bits Eqs. I|14|) are equal to 
— 1 iff an error has occurred on the edge u = (e, e, 1) or 
u = (e, e, d) of the lattice T e . Since each of these errors 
shows itself in a corresponding syndrome bit which is not 
affected by any other error, we can reliably identify these 
errors. This is equivalent to actively correcting them with 
unit success probability. We can therefore assume in the 



subsequent analysis that no errors occur on the edges 
(e, e, 1) and (e,e,d), which concludes the discussion of 
the issue (ii). 

As in [lj| , we define an error chain £ as a collection of 
edges where an elementary error has occurred. Each of 
the two lattices T e and T a has its own error chain. An er- 
ror chain £ shows a syndrome only at its boundary d(£), 
and errors with the same boundary thus have the same 
syndrome. One may identify an error £ only modulo a 
cycle D, £' = £ + D, with d(D) = 0. 

There are homologically trivial and nontrivial cycles. 
A cycle D is trivial if it is a closed loop in T Q (T e ), and 
homologically nontrivial if it stretches from one rough 
face in T Q (T e ) to another. A rough face here is the 2D 
analogue of a rough edge on a planar code [l^. The 
rough faces of T a are on the upper and lower side of C, 
and the rough faces of T e are on the front and back of C 
(recall that no errors occur on the left and right rough 
faces of T e ) . 

Let us now study the effect of error cycles on the iden- 
tification of the state \iP)lr from the measurement out- 
comes. We only discuss the error chains on T a here, which 
potentially affect the eigenvalue Eq. 10. The discus- 
sion of the error chains in T e — which disturb the ZlZr- 
correlations — is analogous. An individual qubit error on 
v G C will modify the X^Xr correlation of lr if it ei- 
ther affects X l , Xr or Xxx ■ That happens if v G Tj^ x ' . 

Now, the vertices in Tj^x correspond to edges in T Q . If 
an error cycle D in T a is homologically trivial, it inter- 
sects in an even number of vertices; see Fig.|3J This 
has no effect on the eigenvalue Eq. ©. However, if the 
cycle is homologically nontrivial, i.e. if it stretches be- 
tween the upper and lower face of C, then it intersects 
Txx m an °dd number of vertices. This does modify 
the eigenvalue Eq. @ by a sign factor of (—1) on the 
l.h.s., which leads to a logical error. Therefore, for large 
system size, we require the probability of misinterpreting 
the syndrome by a nontrivial cycle to be negligible ^3] : 

^prob(f) prob(£ + D\ £ ) « 0. (15) 

£ D nontrivial 

We have now traced back the problem of reconstructing 
an encoded Bell pair \iP)lr. to the same setting that was 
found in 0| to describe fault-tolerant data storage with 
the toric code. Via the measurement pattern Eq. © , we 
may introduce two lattices T a , T e such that 1) Syndrome 
bits arc located on the vertices of these lattices, 2) In- 
dependent errors live on the edges and show a syndrome 
on their boundary, 3) Only the homologically nontrivial 
cycles give rise to a logical error. This error model can be 
mapped onto a random plaquette Z^-gauge field theory in 
3 dimensions [T2I Il5j which undergoes a phase transition 
between an ordered low temperature and a disordered 
high temperature phase. In the limit of l,d — > 00, full 
error-correction is possible in the low temperature phase. 

In our setting, the error probabilities for all edges are 
equal to p. For this case the critical error probability has 
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been computed numerically in a lattice simulation [lfij . 
p c = 0.033 ± 0.001. This value corresponds, via Eq. JSJ, 
toT c = (0.296 ± 0.003) A. 

Remarks: 1) The error model equivalent to Eq. J^J, 
i.e. Z-errors only, is very restricted. We have a physical 
motivation for this model, but we would like to point out 
that the very strong assumptions we have made about 
the noise are not crucial to our result of the threshold 
error rate being non-zero. One may, for example, gen- 
eralize the error model from a dephasing channel to a 
depolarizing channel, with p x — p y = p, = p' '/3. Then, 
two changes need to be addressed, those in the bulk and 
those on the faces L and R. Concerning the faces, the 
errors on the rough faces to the left and right of T e can 
no longer be unambiguously identified by measurements 
of the code stabilizer (|14() . which raises the question of 
whether — for depolarizing errors — it may be these sur- 
face errors that set the threshold for long-range entan- 
glement. This is not the case. To see this, note that two 
slices of 2D cluster states may be attached to the left 
and right of C, at 7/3 = 0, —1 and U3 = d+ 1, d + 2. The 
required operations are assumed to be perfect. They do 
not change the localizable entanglement between the left 
and right side of the cluster C because they act locally on 
the slices — 1..1 and d..d + 2, respectively. The subsets A 
and B of spins are re-located to the slices —1 and d + 2, 
with the corresponding changes in the measurement pat- 
tern. The effect of this procedure is that the leftmost and 
rightmost slice of the enlarged cluster are error-free , 
and only the bulk errors matter. 

Concerning the bulk, note that the cluster qubits mea- 
sured in Z-basis serve no purpose and may be left out 
from the beginning. Then, the considered lattice for the 
initial cluster state has a bec symmetry and double spac- 
ing. The lattices T a , T e remain unchanged. Further, X- 
errors are absorbed in the X-measurements and K-errors 
act like Z-errors, such that we still map to the original 
Z2 gauge model [l^ at the Nishimori line. The threshold 
for local depolarizing channels applied to this configura- 
tion is thus p' c = 3/2p c = 4.9%. In addition, numerical 
simulations performed for the initial simple cubic cluster 
and depolarizing channel yield an estimate of the critical 
error probability of p" c = 1.4%. 

2) Finite size effects. We carried out numerical simula- 
tions of error correction on an I x I x d lattice with periodic 
boundary conditions (as opposed to the open boundary 
conditions of the planar codes within the cluster state). 
For differing error rates below the threshold value of 2.9% 
[Ti| , we found good agreement for the fidelity F between 
the perfect and the error-corrected encoded Bell state 
with the model F ~ exp(— dk\ exp(— Ikz)). Some data is 
shown in Fig. ^corresponding to a Z-error rate of 1.0%. 
Provided that planar codes and toric codes have similar 
behavior away from threshold, our simulations suggest 
that, in order to achieve constant fidelity, the length / 
specifying the surface code need only scale logarithmi- 
cally with the distance d. 

3) For even d, the construction presented above can be 




FIG. 4: This figure plots data for simulations of error correc- 
tion on an L x L x dtoric lattice, with periodic boundary con- 
ditions in the first two directions, for various L and dtoric (d = 
2dtoric + 1,1 — 2L). The error rate is p = 0.01. The logs are 
base e. Two standard deviations above and below the com- 
puted values (as given by statistical noise due to the sample 
sizes) are shown by the error bars. The solid lines each have 
slope one, and they are spaced equally apart. This lends good 
support to the model of fidelity F ~ exp(— dki exp(— Zfe)) for 
error rates below threshold. 



used to mediate an encoded conditional Z-gate on distant 
encoded qubits located on slices 1 and d. 



III. UPPER BOUND 

In this section we analyze the high-tcmpcrature behav- 
ior of thermal cluster states and find an upper bound on 
the critical temperature T c . Our analysis is based on the 
isomorphism between cluster states and the so-called Va- 
lence Bond Solids (VBS) pointed out by Verstraete and 
Cirac in 01 which can easily be generalized to a finite 
temperature. 




FIG. 5: (color online) (a) Correspondence between physical 
and virtual qubits. Domains are shown by dashed lines, (b) 
A bipartite cut of a cubic lattice. The regions A and B are 
highlighted. 

With each physical qubit u G C we associate a domain 
u.* of d(u) virtual qubits, where d(u) = |neigh(w)| is 
the number of nearest neighbors of u (see Fig. (a)). 
Let us label virtual qubits from a domain u.* as u.v, 
v 6 neigh(u). Denote E to be the set of edges of the 



G 



lattice C and define a thermal VBS state p VBS as 
Pvbs = Yl 4 ( 7 + UJeXu 

(16) 

Here {uj e } are arbitrary weights such that < oj e < 1. 
It should be emphasized that p VBS is a state of virtual 
qubits rather than physical ones. Our goal is to convert 
p VBS into p cs by local transformations mapping a domain 
u.* into a single qubit u G C. The following theorem 
is a straightforward generalization of the Verstraete and 
Cirac construction (here we put A/2 = 1). 

Theorem 1. Let p cs be a thermal cluster state on the 
3D cubic lattice C at a temperature T = j3 . Consider 
a thermal VBS state p VB s as * n Eq. l^lb]) such that the 
weights uj e satisfy 

^(it,!)) > tanh (/3) for each u G C. (17) 

u£ncigh (u) 

Then p VBS can be converted into p C3 by applying a com- 
pletely positive transformation VV« to each domain u.*, 

Pes = W(Pvbs), W = (g) W u . (18) 
uec 

Let us first discuss the consequences of this theorem. 
Note that each edge e G E of p VBS carries a two-qubit 
state 

p e = -(I + u e X 1 Z 2 )(I + u> e Z 1 X 2 ). (19) 

The Peres-Horodecki partial transpose criterion 0, 0| 
tells us that p e is separable if and only if uj e < \/2 — 1. 
Consider a bipartite cut of the lattice by a hypcrplane of 
codimension 1 (see Fig. 03(b)). We can satisfy Eq. l(T7|) 
by setting u> e = tanh {(3) for all edges crossing the cut 
and setting ui e = 1 for all other edges. Clearly, the state 
p VBS is bi-separable whenever tanh (/?) < \[2 — 1. But 
bi-separability of p VBS implies bi-separability of p cs . We 
conclude that the localizable entanglement between the 
regions A and B is zero whenever tanh (fl) < \/2 — 1, 
which yields the upper bound on T c presented earlier. 

Remarks: We can also satisfy Eq. I|17fl by setting 
Lu e = uj for all e e E, with uj 6 = tanh(/3). This 
choice demonstrates that p cs is completely separable for 
tanh (p) < (y/2- 1) 6 (that is T « 200). It reproduces the 
upper bound poj of Diir and Briegel on the separability 
threshold error rate for cluster states. 

In the remainder of this section we prove Theorem ^ 
Consider an algebra A u of operators acting on some par- 
ticular domain u.*. It is generated by the Pauli operators 
Z u _ v and X u . v with v G neigh (u). The transformation 



W u maps A u into the one-qubit algebra generated by 
the Pauli operators Z u and X u . First, we choose 

W„(t?) = WlvWu, W u = |0W)(0| + |l® d ("))(l|. 
One can easily check that 

W ] u Z u . v = Z u Wl and Z U . V W U = W U Z U , (20) 

for any v G neigh (u). As for commutation relations be- 
tween W u and X u _ v one has 

(iLjeneigh(u) Xu.v J W u = X u , ^-n 

w u (Tl veS x u . v ) w u = o, 

for any non-empty proper subset S C neigh (u). Taking 
W = <8> ue c ^ u an d usm S Eqs. (J20J, lt2T|l one can easily 
get 

w{pvbs) = n ( j + 7 ?" a '«) > Vu= n w( U)t)) . 

wGC vGneigh (it) 

(22) 

We can regard the state in Eq. l]_!2l) as a thermal cluster- 
state with a local temperature tanh (/3 U ) = r\ u depending 
upon u. The inequality of Eq. (|17fl implies that (5 U > /3 
for all u. To achieve a uniform temperature distribution 
/?« = one can intentionally apply local Z-errors with 
properly chosen probabilities. 



IV. CONCLUSION 

Thermal cluster states in three dimensions exhibit a 
transition from infinite to finite entanglement length at 
a non-zero transition temperature T c . We have given a 
lower and an upper bound to T Ci 0.3 A < T c < 1.15 A 
(A = energy gap of the Hamiltonian) . The reason for T c 
being non-zero is an intrinsic error-correction capability 
of 3D cluster states. We have devised an explicit mea- 
surement pattern that establishes a connection between 
cluster states and surface codes. Using this, we have de- 
scribed how to create a Bell state of far separated encoded 
qubits in the low-temperature regime T < 0.3 A, making 
the entanglement contained in the initial thermal state 
accessible for quantum communication and computation. 
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